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Abstract 

With tangent angle perturbation approach the axial symmetry deformation 
of a spherical vesicle in large under the pressure changes is studied by the 
elasticity theory of Helfrich spontaneous curvature model. Three main results 
in axial symmetry shape: biconcave shape, peanut shape, and one type of 
myelin are obtained. These axial symmetry morphology deformations are 
in agreement with those observed in lipsome experiments by dark-field light 
microscopy [Hotani, J. Mol. Biol. 178, 113 (1984)] and in the red blood cell 
with two thin filaments (myelin) observed in living state (see, Bessis, Living 
Blood Cells and Their Ultrastructure, Springer- Verlag, 1973). Furthermore, 
the biconcave shape and peanut shape can be simulated with the help of a 
powerful software, Surface Evolver [Brakke, Exp. Math. 1, 141 (1992)], in 
which the spontaneous curvature can be easy taken into account. 

PACS numbers: 87.16.Dg, 47.10.+g, 68.15.+e, 02.40.Hw 
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I. INTRODUCTION 



Since the discovery of optical microscope, people have recognized that the human red 
blood cells (RBCs) QTJ] and lipsomes can be in many beautiful and strange shapes depend- 
ing on the environment in which they exist, such as pH value, osmotic pressure difference. 
The simplest model to explain these shapes is viewing RBC and lipsome as a vesicle of 
amphiphile bilayers. 

For a long time, both the experimental and the theoretical aspects of amphiphile bilay- 
ers and monolayers have gathered much attention from physicists and chemists ||. Helfrich 
has shown the importance of spontaneous curvature and has developed a phenomenological 
theory for the elasticity of fluid membranes by an analogy with the curvature elasticity of 
liquid crystals M. Many authors have studied the deformation of vesicles with different 
shapes among which the spherical vesicle plays a special role in these studies. The general 
shape equation of vesicles has been derived with differential geometry method ||, then some 
interesting numerical and analytical results have been reported |||7|]. In the equation, the 
osmotic pressure difference between the external and internal sides of the vesicle membrane, 
especially the spontaneous curvature, and elastic modulus are important quantities. Among 
them, the influence of pressure on morphologic appearance of vesicle is the most important 
one to be considered. In shape study, RBC becomes an attractive model of vesicle in living 
cells. On the other hand, Hotani reported the lipsome shape transformation pathway when 
it was put in solution at certain concentration, and pointed out that osmotic pressure was 
found to be the driving force for the sequential transformations @. In the experiment the 
concentration changes slowly, and so does pressure, therefore, the shape formation of the 
lipsome can be seen as an equilibrium problem. In Hotani's observation, single sphere finally 
formed and seemed to be stable forms because they showed no further morphological change. 
Jenkins and Peterson studied in detail the problem of the stability of RBC shapes ||, in 
which spherical vesicle was limit case. Ou-Yang and Helfrich studied instability and defor- 
mation of a spherical vesicle by pressure, and they found that any infinitesimal deformation 
corresponding to spherical harmonics Y\ m would require a pressure difference larger than 
some threshold values ]9[]. The nonaxisymmetric deformation of RBC is investigated with a 
contour perturbation approach [1(J. Obviously, the above approach is not able to deal with 
a large deformation of a vesicle shape , such as the myelin form in RBCs [fjD]. It then becomes 
a challenge to search for a method by which the large deformation of a spherical vesicle can 
be treated in theory of perturbation. In the present work we report an useful issue on solving 
the problem. The axial symmetry vesicle contour is expressed as Z = — J tan ^>dp where \& 
is the tangent angle of the contour (see Fig. [I]). One can find that a small change of $ at 
^ — > 7i / 2 is able to cause a large deformation in contour. Just like the method of Deuling 
and Helfrich ||, we introduce an alternative variable of sin\l/, u, thus, the contour Z will 
change with considerable amplitude as sin \1/ has small perturbation, with the issue we can 
calculate the large deformation of the contour via perturbation approach. 

In this paper, we do some calculations about the large deformation of a spherical vesicle 
with the new approach via pressure perturbation, varied biconcave shapes and one kind of 
peanut shape can be obtained. Especially we show one type of myelin form with two long 
filaments attaching on the spherical vesicle. To confirm the analytic calculation we also 
perform some computer simulations with the help of a powerful software, Surface Evolver 
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II. PERTURBATION THEORY 

Equilibrium shapes of phospholipid vesicles are assumed to correspond to the local min- 
ima of the elastic free energy of these systems. In the Helfrich spontaneous curvature model, 
the spontaneous curvature plays a fundamental role in accounting for the different morpho- 
logic appearances of vesicles ||]. According to Helfrich's theory, the free energy of a vesicle 
is written as 



F = ht c j{c x + c 2 - c ) 2 dA + Ap J dV + A j> dA, 



(1) 



where k c is the bending elastic modulus; c\ and c 2 are the two principle curvatures of the 
surface of the vesicle, cq is the spontaneous curvature to describe the possible asymmetry of 
the outer and inner layers of membrane, and dA and dV are the surface area and the volume 
elements for the vesicle, respectively. The Lagrange Multipliers Ap and A take account of 
the constraints of constant volume and area, respectively, and can be physically understood 
as the pressure difference between the external and the internal environment and the surface 
tensile coefficient. The general shape equation has been derived via variational calculus || 
to be 



Ap - 2XH + k c (2H + c )(2H 2 - 2K - c H) + 2k c y 2 H 
= 0, 



(2) 



where H and K are the mean and the Gaussian curvatures, respectively, and y 2 is the 
Laplace-Beltrami operator ||. Assuming that the shape has axial symmetry, the general 
shape equation becomes a third order nonlinear differential equation 
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where p is the distance from the symmetric Z axis of rotation, ^(p) is the angle made by the 
rotational axis and the surface normal of the vesicle, also the tangent angle of the contour 
Z(p) (see Fig. |1|). Let sin\l/ = u, p = x the equation takes the following form 
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where prime means derivative of u(x) with respect to x. It is obvious that a sphere with 
radius r is always a solution of Eq. @) (it = x/r ), and its pressure difference = Ap 
must obey 



2\r 2 Q + fc c c r (2 - c r ) = 0. 



(5) 



Jenkins [f|, Ou-Yang and Helfrich || have studied the stability of the spherical vesicles. 
Now we calculate the first order pressure perturbation contribution to u using the Eq. (f|). 
When the pressure difference Ap is slightly deviated from its equilibrium value Apo, i.e. , 
Ap = Apo + SAp, Eq. (|5]) is no longer satisfied, what will be the shape of the vesicle? Here 
we give an answer to this question with our detailed calculation. 

By expanding u as u = u + U\ +u 2 + ... with u\ ~ SAp, u 2 ~ (5Ap) 2 , ... , while u takes 
x/tq, we find u\ to satisfy the following linear third-order differential equation 




5Ap 

k r 



0. 



(6) 



where 5Ap is the pressure perturbation. After introducing a new variable £ = x/tq, i.e., uq 
takes £, the above perturbation equation becomes 



:i -ex 



+ &x 



+ (^3 -7)«1+P=0, 



(7) 



where two dimensionless parameters are: 6 = r^X/kc + CqtI/2 — 2c r , p = r^Ap/k^ and 
now prime means derivative of u(£) with respect to £. It is obvious that u\ = (p/2b)£ is a 
particular solution. Here we suppose that 6^0. For the case of b = 0, we discuss it in the 
next section. Then the first step for solving Eq. (|7|) is to solve the following homogeneous 
equation 
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It is lucky that u% — £ 1 is a particular solution of this homogeneous equation. We denote 
it as itio = £~\ and let 
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ui = - I z(£)dS, (9) 
Eq. (||) is thus reduced into the following second-order homogeneous differential equation 

_ e v - & + a - he)z = o. (fo) 



To solve it, we refer to the following standard form of the equation fll4 

x 2 (ax n - l)y" xx + x(atx n + q)y' x + (arx n + s)y = 0. 

There are four roots Ax, A 2 and B 1 , B 2 of the quadratic equations: A 2 — (q + l)A — s = and 
B 2 — (t — 1)B + r = 0. Then four parameters c, a, (3, and 7 can be defined by the relations: 
c = Ai, a = Al ^ Bl , P = 7 — Al ^ A2 + 1 an d the solution of the standard equation has 

the form y = x c u(ax n ), where ct>(£) is the general solution of the hypergeometric equation 

£(£ - + [(a + /3 + 1) - 7]^ + a(3uj = 0. 

In the case of Eq. ([TOD , we have a = 1, n = 2, t — 0, q = 1, r = b and s = —1. Two 
quadratic algebra equations associated with Eq. ( |10|) are introduced: A 2 — 2A + 1 = and 
B 2 + B + b = 0. Their roots are: A 1 = A 2 = 1, and fii, B 2 = (-1 ± VI - 46) /2. From them 
three characterized numbers can be calculated : 
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Ax + Bx 1 + vT^46 
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7 = 1. (11) 



Then we have 



*(0 = &>(e a ), (12) 

where a;(r/), 77 = £ 2 , is the solution of Gauss hypergeometric equation: 

r](r] - IV" + [(« + /? + 1)77 - 7 ]u/ + «/3cj = 0. (13) 
At once, let its solution uix take the hypergeometric function WA 



as 



u 1 = F(a,P,<y,7 t ). (14) 
Because 7 = 1 is an integral, we have the second solution of Eq. (|T3| ) |15 | 

^2^) = F(a,(3,7,r]) ln(r^) 

+ E (fe!)2 V {</>(<* + k) 

+4>(f3 + k) -20(1 + A;)}, (15) 

where (A) = 1, (A)„ = A(A + 1)...(A + n - 1) = (n > 1), 0(r) = ^» , and T(r) 

is Euler Gamma function. Finally, we find two solutions of Eq. ( |I0| ) 
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z x = £0^(77) = ^i(e 2 ) = £F(a, 0, 1, e), 
Z2 = & 2 (r)) = ^ 2 (£ 2 ) = f F(a, P, 1, £ 2 ) ln^ 2 + 
g(^*^ (o + fc) 

+ fc) -20(1 + A;)}, (16) 

and two solutions of Eq. @ m u = tui(£, 2 )d£, 2 , w 12 = ^ f cu 2 (£, 2 )d£, 2 . To obtain the 
detailed forms of ttu and U\ 2 two integrals, / F(a, P, 7, 77)^77 and / F(a,ft, 7,77) In 77^77, need 
to calculate. We have (see the details in appendix) 
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1 
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where F'(a, P, 1, 77) = -j-F(a, P, 1, 77) = aPF(a + 1, /3 + 1, 2, 77). Now, we can write down 
three independent solutions of Eq. (§): 
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+ (a + /3-l)£ 2 F(a,/3,l,a ln£ 2 + 

2-a-p 



(a + (3 - a(3 - l) 2 
[£ 2 F(a,/3,U 2 )(a + /3-l) 
+£ 2 (£ 2 -l)a/3F(a + l,/? + l, 2, £ 2 )]} 

fc=0 

-20(1 + k)}. (19) 

From Eq. (^), Eq. (Q) and Eq. (|i~9"|) we get the general solution of Eq. (||) as 

1 1 

M = M + PCIQ- + pen Un + PC12W12 + P^£, (20) 

where cio, en, and C12 are three integral constants. 

In general, the hypergeometric function can be expressed as the series 

F(a,/5, 7 ,,) = l + g^--, (21) 

where (A) = 1, (A) n = A(A + l)...(A + n — 1). F(a, /?, 7, 77) is a fortiori convergent for \rj\ < 1. 
In the present case 7 = 1, a = (1 + vl — 46) /4 and /3 = (1 — yl — 46)/4, the condition 
_Re( 7 — a — /3) = | > is satisfied, so F(a,/3, l,rj) is conditional convergent while 77 = 1, and 
F(a,/3, 7, 1) = fr~^^E^] • But this case is not true to F(at + 1,(3+ 1, 2, 77) when rj = 1, 
and the analytical continuation for it gives the complex result. So we limit our discussion 
within < £ < 1 to avoid the divergence and imaginary number. From another point of 
view, if we consider that £ can be or be more than unit, that three integral constants cio, 
en and C12 in Eq. fllTf ) must be zero, so the solution of Eq. (f|) will be uo (sphere solution) 
adding that particular one, i.e., 

P f 

the sphere still keeps its shape under a pressure perturbation varying its size only. After 
obtaining the total solution of Eq. (§), using the integration 

tan^dp, (22) 

we can give the contour of the axial symmetry shape of vesicle with deformation. Our 
numerical calculation about Gauss hypergeometric function is worked out with the famous 
software, Mathematica. In the following sections some cases are shown for details. 



III. SPECIAL CASE 

First we discuss the special case of 6 = in the Eq. ([/]). That is to solve the following 
equation 
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(1 - £ >"' + (| - 3£K - i«' + iu + p = 0. (23) 

It is convenient to check that u = p/2£ log£ is a particular solution. So the task remained 
is to solve the associated homogeneous third order equation 

(1 - e)u'" + (| - 3£K - j 2 u' + ±u = 0. (24) 

It is easy to verify that U\ — £, u 2 — l/£ are two independent solutions. According to the 
handbook [14j], the general solution of Eq. (^4j) is written as 



u = c\Uiq + c 2 u 2 + c 3 (m 2 y U\9dx — u\ J u 2 6dx) (25) 

where 

9 = exp[ / {h/ h)dx){uxu' 2 - ^ns) -2 , 



The final result is 



/ 2 = 2/e-3e, ^ / 3 = 1 - 



u = i + Cif + — + c 3 (£ log ) 

+f£log£- (26) 



If c 2 = c 3 = 0, we have 



« = e+cie+feioge (27) 



Eq. (|27|) satisfies the general shape equation Eq. (|[). It can give the biconcave shape @ and 
many interesting shapes described as have been shown in |To] . 



IV. BICONCAVE AND PEANUT SHAPE 

With the help of microscopes people have made detailed observation of RBC which can 
assume various shapes 0. Generally, RBC takes biconcave disk shape (in blood capillary it 
has very large deformation [ITU), while some pathologic cells take other abnormality (such 



as in sickle cell disease). If RBC is subjected to different environment, various pH values 
for example, it will make some considerable deformations, such as, Cup (Stomatocyte), Bell 
(Codocyte), Sea urchin (Echinocyte) et ai. Moreover, beginning with biconcave lipsome 
vesicle Hotani present many beautiful transformations among which the peanut shape was 
included ||. Surely the biconcave shape attracts many attention. In this section we use the 
general theory of the section II to show that spherical vesicle can transform into biconcave 
shape via pressure perturbation. Because ^ = at p — and 1 I' = 7r/2atp = p m (see 
Fig. D), there are two boundary conditions at p = and p = p m , i.e., £ = and £ = £ m 
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u|£=o = 0; u\t =u = 1, (28) 

where £ m = p m /r . Thus there are two relations between the three coefficients c w , Cn and 
Ci 2 as 



ClO 



-Cl2 



2(a + /3 - a(3 - 1) 



2-a-p 1 
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r( 7 )r( 7 -g-/3) i , i , 
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r( 7 )r( 7 -«-/5)(« + /3-i) ' 



(29) 



So far, in the above solution there is still one relation to be determined. We can calculate it 
via using the conservation of area of sphere. For a rotational surface the area element can 
be expressed as 

ds = 2vrx(l - u 2 )~ l/2 dx. 

Using reduced dimensionless area s r defined by ds = 47ir 2 ds r , we have ds r = |^(1— m 2 )^ 1 / 2 ^, 
and another boundary condition to take account of the area conservation 

' £(1 - u 2 )-^ 2 dC = I. (30) 

Basing on Eq. (|20"D , Eq. (|29| ) and Eq. (|50]), the values of cio, cn and c\2 can be fixed 
completely. So the shape is also determined with Eq. (^2|). But in practical computation the 
value of £ m in Eq. ( |3"0"D can not be directly predicted. Therefore, in our computation, we treat 
it as a free boundary problem, i.e., we let £ m takes a series of value, correspondingly, a series 
of deformation of the original spherical vesicle and their total energy can be calculated ||TH[ . 
The less total energy, the easier the deformation can be observed in practice (experiments). 
From the Eq. ( pU|) the different shapes can be got via Eq. (P2j). We show our numerical 
results in Fig. 0, Fig. |3| and Fig. [|. In fact, in Fig. |2| only the top right quadrant contours 
of vesicle are shown. In order to obtain the total vesicle surface two operations need to 
be done in succession: at first, rotate it with the Z axis, this give the upper part of the 
surface; second, take mirror reflection. In order to describe the parameters corresponding 
to the curves in Figs. ^|-6, we define a new parameter v as v = V/Vq where V is the volume 
of vesicle and Vq is the volume of original spherical vesicle. The pressure corresponding to 
each curve is defined as Ap ~^ po , Ap is the pressure related to the deformed vesicle while 
Ap is the pressure difference of the original spherical vesicle. The shapes given in Fig. |2| 
are obtained at the constant spontaneous curvature and surface tension, i.e., cqTq =-2.4 and 
\r 2 /k c = 0.01. The values of ( Ap ~ Ap ° , v) are calculated as (-lO^ 6 , 0.99), (0.01, 0.3), (0.01, 
0.5), (0.05, 0.3), (0.01, 0.62) for curves 2-6, respectively. The curve 1 is the initial shape 
of sphere. Two three-dimensional shapes of biconcave and peanut are shown in Fig. |3| and 
Fig. ||, respectively. In principle, energy corresponding to every curve can be calculated with 
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the formula in |1| . The above calculated biconcave shape can be simulated on the basis of 



the elasticity theory of Helfrich spontaneous curvature model by Surface Evolver ||11|| . The 
simulated results, corresponding to biconcave and peanut shape, are shown in Fig. [5] and 
Fig. |^, respectively, in which the A/A gives the ratio of surface area of vesicle to that of 
the initial spherical vesicle. Both they are observed in blood cell |1| and lipsome vesicles [Q. 
Furthermore, We turn to analysis the behavior of deformation near rotational axis. After 
obtaining the general solution Eq. fl2"U|), the behavior of deformation near the rotational axis 
can be got by expanding it to first order of p as 

1 1 -ad r> 
u = £ + -pc^ + pc 12 (ln£ + - + 2 — R - )i + j-X 

2 2 a + p — ap — 1 2b 

= A£ + Bfln£. (31) 

This equation describes the shape near polar points, and the equation is very similar to the 
biconcave shape of red blood cell JF|. In this equation there is no any singular point in the 
region < £ < 1. 



V. MYELIN FORM 

Myelin form may originate from all blood cells Jl|. Shadowing technique in microscopy 
reveals that myelin forms are hollow structures. When a RBC is in aging, aged and damaged 
states, it gives rise to large myelin forms which may take various types: filaments, beads 
or strings of beads. These filaments, which are easily seen with phase contrast microscopy, 
may remain attached to the surface of the RBC at one end. Here we choose the so called 
Delaunay's surface solution |7] from the general solution of the perturbation Eq. (||), which 
can give a type of myelin form generating from the vesicle. In the book of Bessis pi, there is 



a photo of this kind of Medusa head. Now let us first discuss the general solution Eq. ( |20| ) 
in the case of en = Ci 2 = 0: 

u=(l + ^)e + <V£. (32) 

Substituting it into the general shape equation Eq. ([|) yields two relations: l+p/2b = cor /2, 
and Ap = cq\. It should be noticed that Ap = Ap + p. Ap Q can be solved from Eq. (^j), 
then we find 

p = (2c - clr )r -^1 + 2^1. ( 33 ) 

for ru/- 



Considering the mentioned relation 1 +p/2b = CqTq we have 

(clr - 2c )r + ^° - ^1 = 2 b - c r b. 



The value of cio in Eq. Q3"2"D can be determined via boundary condition: u\^ m — 1 (£ = £ 



in ■ 



^ = 7r/2). It leads to c w = (£ m — (1 + ^)£m) an( i finally we obtained from Eq. fl32|) the 
Delaunay's solution 
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u= c^ + U-c^0J2_ (34) 

With different cor and £ m , the various contours can be numerically calculated using Eq. (|22|) 
(see Fig. [7] and Fig. ||). Not total length of filament is drawn in this two figures. In fact, those 
filaments are swollen distally and our present calculation can not give an exact description 
yet. 



VI. CONCLUSIONS AND DISCUSSIONS 



Now we give here our main conclusions as what follows: With tangent angle perturba- 
tion approach, we can calculate the large deformation of spherical vesicle under pressure 
perturbation. From the general perturbation solutions, the biconcave and peanut shapes 
can be obtained and a kind of myelin form is shown to be existed. Our all calculations are 
based on the elasticity theory of Helfrich spontaneous curvature model and it does give good 
accordance with some complex shapes of lipsome vesicles (see photographs in f2|) and RBCs 
(see photographs in []]]) with computer simulation |T§|,|I|,|2| . 
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APPENDIX: TWO INTEGRALS RELATING TO GAUSS HYPERGEOMETRIC 

FUNCTION 



In the appendix we give the derivation of two relating integral about Gauss hypergeo- 
metric function. The main idea and method are to use hypergeometric equation itself. From 
the hypergeometric equation Eq. ( |T3| ) 



we have 



and 



x (1 - x)u" + [7 - (a + (3 + l)x) J - a/3u = 0, 



u = \x (1 - x)u" + [7 - (a + 13 + 1)] u/, 

ap 



x(l-x)u' 7-I+ (a + P + l)x /A1 . 
udx = ^ - 1 + — — ^ - ' u. Al 

l-a-(3 + a(3 I - a - (3 + a(3 v ' 

Using the same method and with the key point 7 = 1, through a lengthy derivation, we 
obtain 



11 



J In xoudx = 



a + f3 + 1 
[uj'x (x — 1) lnx + (a + P — l)xu\nx 
/-, \ i 2 — a — (3 

(a + (3 — ap — 1) 
[iu (a + (3 - 1) + x(x - l)oj']. 

(A2) 
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FIGURES 




FIG. 1. Convection of angle Vl/ for the contour of axial symmetry vesicle. At the north pole the 
^ value takes zero. p m gives the boundary position of p. 




FIG. 2. The numerical results from the solution of Eq. (|20[). Only the top right quadrant 
contours are shown. Curve 1 corresponds to spherical vesicle, while curves 2-6 show much swelling 
or oblate biconcave shapes. The energy of each vesicle can be calculated with the formula given in 
[18]. 




FIG. 3. Biconcave shape in three-dimensional (half of upper part is shown) with £ m = 0.9999, 
v = 0.1 (the ratio of volume of the deformed vesicle to that of the initial sphere), Xr^/kc = 0.01, 
co r o = —3.0. When the parameters are adjusted, the other biconcave shapes can be formed, much 
oblate or much swelling (see, Fig. [2]). 
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FIG. 4. Numerical result of upper-half peanut shape in three-dimensional with £ m = 0.8, 
v = 0.4, Xrl/k c = 0.01, and c r = 0.01. 




FIG. 5. Biconcave shape in three-dimensional simulated with Surface Evolver. The parameters 
take Hq=-1.38, k c =l, A/Aq=0.96 (the ratio of surface area of vesicle and that of the initial sphere), 
v=0.92, and Arg//c c =0.48. 




FIG. 6. Peanut shape in three-dimensional simulated with Surface Evolver. The using param- 
eters are co?"o=1.38, k c =l, ^4/Ao=0.83 (the ratio between surface area of vesicle and that of the 
initial sphere), v=0.79, and ArQ//c c =0.48. 
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FIG. 7. The numerical results corresponding to Eq. Only the top right quadrant contours 

are shown. The parameters corresponds to curve 4 (contour in Fig. ||) are: cqTq = 2, £ m = 0.9. 




FIG. 8. The calculated shape of upper-half of the myelin form of spherical vesicle in 
three-dimensional with c^vq = 2, and £ TO = 0.9. 
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